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| 7j I Abstract 

■ The effect of residual interaction between nucleons (quasiparticles) on shell oscil- 

, lations of the masses of spherical nuclei is considered. The singularity of the ground 

^ I state energy of the system in the vicinity of nucleon magic numbers is analyzed for 

' various types of the dependence of residual interaction on orbital momentum of the 

. quasiparticle. It is shown that only the perturbation band width of the Fermi distri- 

^—^ \ bution due to residual interaction which is proportional to the square of the angular 

' momentum vector is consistent with the character of the magic cusps. The coupling 

. constants between the quasiparticles are determined on the basis of the available data. 

I The constant decreases rapidly with increase of nuclear radius. Possible consequences 

^ ' pertaining to the energy spectrum of infinite nuclear matter are discussed. 

^ '. 

^ : 

Ch ■ 1 Introduction 

I In a preceding paper [1] we considered the influence of the sharpness of the boundary of 
I quasiparticle Fermi distribution on the characteristics of spherical nuclei. It turned out that 
■ - - ' the number N{pf) of single-quasiparticle states located below the boundary pf = kfR ^ 1 
{kf is the corresponding value of the wave number and R is the radius of the nucleus) is 
strictly speaking not an analytic function of its argument. In the p scale, the singularities 
of the function are equidistant, with an interval 7r/2; a qualitative idea of their character 
is gained from Fig. la. These are the magic cusp points (jumps of the first derivative); a 
consistent analysis shows that singularities of this type are possessed also by the energy of 
the ground state of the body, i.e., the mass of the nucleus. Figure 2 shows schematically a 
typical experimental plot of the mass in the vicinity of a magic nucleus ^. We have in mind 

*Zh. Eksp. Teor. Fiz. 63, 1961-1977 (1972) [Sov. Pys. JETP 36, No. 6, 1036-1044 (1973)] 
"'^We emphasize the physically unrealistic character of the "sub-magic" nucleon numbers, which can be 
formally set in correspondence with the filling of each j-level in a certain spherically-symmetrical potential 
well. No phenomenon pertaining to the nucleus as a whole is observed experimentally in this case; nor does 
a consistent theoretical analysis seem to predict any phenomenon. Moreover, not all the characteristics of 
the referred-to "well" admit of a clear-cut physical definition. The question is: what takes place at the 
point where the magic nucleus is located? Does the Fermi boundary rise because the occupation of the last 
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relatively "weak" singularities; they appear in the higher-order terms of the expansion in the 
reciprocal powers oi pf. We note that in this approximation, which is of interest in nuclear 
physics, there are no grounds whatever for identifying the function N{pf) with the total 
number of true particles N{pf). Furthermore, we can advance the following considerations: 
The roughly intuitive concept of the quasiparticle as moving in the average field produced 
by practically all the particles of the body accounts nevertheless quite well for the main gist 
of the phenomenon. Since the dimensionless parameter Pf = kfR is determined, roughly 
speaking, from the characteristics of this average field, it is natural to assume this field to 
have a smooth dependence on the total number of particles N, as is indeed assumed (see 
Fig. lb). It is easy in practice to make the transition to the physically most significant N 
scale, for example, by expanding in reciprocal powers of pf (see [1], formulas (3), (26), and 
(27) and the explanations pertaining to them). 

We shall now comment on the nature of the shell and magic oscillations in spherical 
nuclei from a point of view that reveals clearly the role of the quantum numbers of the 
individual quasiparticles. Of particular importance is the conservation of the orbital angular 
momentum /; the eigenvalues p = kR can then be graphically represented by points on the 
(/,p) plane. It is easily seen from Fig. 3 how they are grouped in the region of relatively 
small orbital angular momenta. The "Rcgge trajectories" drawn in accordance with the rule 
2n + / = p (n is the principal quantum number and p is the number of the trajectory) lie 
near their maximum. More concretely, their form is given by the equation 

Ap=-(/ + l/2)V2p (1) 

(see also Appendix 1). The vertical distance between the curves is equal to n/2, correspond- 
ing precisely to the aforementioned interval between the neighboring shells in the p scale ^. 
We can now easily visualize the situation that arises if the Fermi boundary 

P-Pf (2) 

is assumed to move, say, upward. After the levels of the last Regge trajectory are exhausted 
at the tangency point (see Fig. 3), the density of states dN /dp decreases by jump. A more 
formal natural interpretation is that the function N{p) has an oscillating component Ni(p), 
the period of which is determined by the distance between the trajectories in Fig. 3. 

This specific grouping of the single-quasiparticle levels, due to the quantum number /, is 
refiected also in the energy E of the nucleus. Its oscillating component is equal to 

EiiPf) = -eN,{pf), (3) 

j-level becomes filled (the most frequently advanced point of view) or, to the contrary, does the edge of the 
"well" drop by a jump in a direction opposite to the Fermi boundary? This question has apparently no 
sufficiently distinct physical meaning, all the more since the position of the bottom of the "well" is not a 
rigorous quantitative concept, owing to the strong damping of the deep quasiparticles. These jumps of the 
chemical potential occur by far not after each filling of the next j-level. There are no sub-magic phenomena 
in spherical nuclei at all, there are only true magic nuclei. Their position is given theoretically by formulas 
(19) and (3) of the preceding paper [1]. 

^In connection with Fig. 3, there is a curious, simple, and experimentally confirmed consequence: each 
nuclear shell contains cither one s or one p state, and the corresponding single-particle levels lie close in 
energy to the end of the filling of this shell. For the atom, there is no such theorem. The absence of any 
far-reaching analogy with the shell structure of the atom was already indicated in the preceding paper (see 
footnote 2 of [1]). Certain features of the nuclear structure, which are weakly pronounced in the case of the 
atom, were mentioned also in footnote 1 above. 
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where Ni(pf) is the oscillating part of the number of the states filled with fermions (quasi- 
particles). The proportionality coefficient —e is the first functional derivative with respect 
to the variation of the distribution function near the Fermi boundary; in other words, it 
is the value assumed here by the energy of one quasiparticle. Since, on the other hand, 
e — —dE/dN is the binding energy of the nucleon in the nucleus, we are thus dealing with 
a quantity that can be directly determined from experiment. For case (2) of an absolutely 
sharp Fermi boundary, the form of the function Ni{pf) was determined in [1]. Of particular 
interest is the decrease Ae in the binding energy of the nucleon which is observed in the 
vicinity of the magic nucleus. In the ideal case (2) we have 

(A.)„ = .-f. (4) 

where s is the arithmetic mean of the values of e on both sides of the magic nucleus. 

To attempt to make clearer the correspondence with the level distribution with respect 
to the momenta, shown in Fig. 3, let us consider the following example, which incidentally 
is somewhat artificial: let the quasiparticles fill only all the vacancies that lie below the 
boundary, the equation for which is 

p{l) =pf- g{l + l/2f/2pf. (5) 

Then calculations similar to those in [1] yield 

As = M' (6) 

Thus, as (7 — > 1, when the curvilinear boundary (5) of the statistical distribution coincides 
with the successive trajectories in Fig. 3 (see also Eq. (1)), the oscillations tend to infinity 
in the considered approximation. 

The example (5) illustrates only the large sensitivity of the oscillations to redistributions 
of the quasiparticles with respect to the quantum states; it has, of course, little in common 
with the real situation in spherical nuclei. Finite nuclear dimensions mean apparently that 
it is impossible to have a canonical transformation to quasiparticles with respect to which 
the case (2) of an ideally abrupt "stepwise" Fermi distribution would be realized. This 
circumstance is frequently called in nuclear physics the "residual interaction" between the 
nucleons. It is felt intuitively that the "smearing" of the Fermi boundary, due to this 
interaction, suppresses, generally speaking, the oscillations. It is seen already from Fig. 3 that 
as the "smeared" Fermi boundary moves forward the Regge trajectories, roughly speaking, 
are gradually captured. Indeed, in the experiment the magic jumps Ae are always smaller 
than the ideal value (Ae)o calculated from formula (4). However, by far not any smearing of 
the Fermi boundary is capable of corresponding, even qualitatively, to the character of the 
experimental data. For example, a statistical distribution of the temperature type 

= exp{(p-p,)/r} + l 
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yields for Ni(pf) an expression that is everywhere differentiable any number of times. Ac- 
cordingly, in this case the function Ni{pf) is analytic and there are no magic singularities ^. 

Let us touch also on the question of the most advantageous way of specifying and investi- 
gating the residual interaction. One can imagine, of course, a situation wherein the residual 
interaction between quasiparticles is specified in the form of a corresponding Hamiltonian; 
one such model example will be calculated in the next section. It is important, however, 
not to lose sight of the following circumstance: although the presence of the interaction 
causes, strictly speaking, the energy of the individual quasiparticle to be no longer a definite 
quantity, this does not influence the applicability of formula (3) in the approximation of 
interest to us. Indeed, let the Fermi-distribution smearing, due to the residual interaction, 
have a width Se; this can be naturally interpreted as an uncertainty, of the same order of 
magnitude, in the quasiparticle energy. On the other hand, in the p scale (see Fig. 3, and 
also [1]), the quantity Sp characterizing the oscillations is of the order of unity. Therefore in 
the region of importance for the oscillations we have 5e ~ {d£/dp)5p ~ ^/p/- Recognizing 
that p/ ^ 1, we have 

5£ < £, (8) 

i.e., the first factor in the right-hand side of (3) remains definite with sufficient degree of 
accuracy. In other words, the redistribution of the quasiparticles over the states, due to the 
residual interaction, should be determined in principle from the condition of the minimum 
of the energy of the nucleus as a whole. But then the oscillations can be calculated without 
taking into account the additional energy of the interaction between the quasiparticles. This 
curious feature of the theory will be illustrated by a concrete example in the next section. 

A "dynamic," so to speak, treatment (meaning that the interaction Hamiltonian between 
the quasiparticles is explicitly specified) would be highly ambiguous. Since, in any case, there 
is no canonical transformation that leads to an absolutely abrupt Fermi boundary (2), we 
apparently have no sufficiently reasonable criterion for a unique choice of the transformation 
to a new quasi-static Hamiltonian. In accordance with (8), the residual interaction can be 
more adequately described by directly specifying the quasiparticle distribution function with 
respect to the states. By regarding this function w(p, /) in a certain sense as a primary con- 
cept, we can apparently hope to use successfully its simple single-parameter approximations. 
Such a characteristic feature of the phenomenon as the tendency of the residual interaction 
to decrease with increasing dimensions of the nucleus is likewise fairly well accounted for in 
this case. We shall return to the pertinent questions in the last two sections. 

2 Simplest model of residual interaction 

When choosing a model example, it is desirable to take into account the experimental fact 
that the spin of an even-even nucleus in the ground state is equal to zero, and that for odd 
nuclei the spin always has the single-particle value within the framework of the shell model 

^The last statement can be regarded also as a consequence of a theorem of more general character. It is 
easy to verify that any coiitiimoiis distribution of quasiparticles over the states, whi(;h docs not depend on /, 
results in no singularities of the type of a jump in the first derivative. In the next two sections of this paper 
we obtain a relation compatible with the experimentally observed picture of the magic phenomena between 
the width of the smearing of the Fermi distribution and the quantum number /. 
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(see, for example, [2]). This suggests the expression 



Hlnt = -Gj a^,at^,a_mam (9) 

m,m'>0 

for the interaction between quasiparticles pertaining to the same j-lcvel. Here and are 
the creation and annihilation operators of a quasiparticle with a 2;-projection of the angular 
momentum equal to m. The Hamiltonian (9) can be diagonalized exactly (see, for example, 
[3]); the eigenvalues are given by the well-known formula of Racah and Mottelson 

E,^, = -G,h,[n,-{h,-l)-s,l (10) 

where 2Vtj = 2j + 1 is the total number of vacancies, hj is the number of interacting pairs at 
the j-level, and Sj is the number of noninteracting quasiparticles (seniority). 

It is clear that the smallest energy of the system corresponds to zero seniority, and 
consequently for the nucleus as a whole the situation reduces to a minimization of the sum 

E = J][2£,-6,. - GM^^ - bj + 1)] (11) 

j 

over the j-levels {ej is the initial value of the quasiparticle energy). The redistribution of 
the quasiparticles should be visualized as occurring at a zero variation of the quantity 

N = Y^2b,. (12) 
j 

This additional condition can be easily taken into account by the method of Lagrange mul- 
tiplier. Taking also the Pauli principle into account (i.e., the condition < 6^ < flj), we get 

f 1, £,-£/<A_, 
_ _ ) Gj{ilj+l)-2{ej-e}) A p ^ A (^'r\ 

- - S 2G~Qi ' < -^f < K^'^) 

^ [ 0, ej-ef>A+, 

where Ej is the chemical potential and Azp = ^Gj{l ^Qj). We present also the interaction 
energy (10) corresponding to the equilibrium distribution: 

— GjQj, Ej — Ef < A_, 

Eint - \ iGj ' < ^3 -^f < 

0, Ej -Ef> A+. 

We note that the initial expression (10) was the result of simultaneous diagonalization of 
the Hamiltonian (9) and the operator of the number of quasiparticles at the j-level; then 
the quantum numbers bj are naturally integers. However, in the subsequent derivation of 
the distribution (13) it becomes necessary to differentiate already with respect to variables 
bj. They are consequently determined with an error of the order of unity. This accuracy 
corresponds precisely to the macroscopic character deduced in [1] for the investigated phe- 
nomena. We ultimately replace flj ± 1 by f2j = j + 1; with the same accuracy, we can replace 
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the total angular momentum j of the quasiparticle by its orbital momentum I and neglect 
the spin-orbit interaction 

We now change to more convenient variables, using a notation close to that of [1] (see 
also [4]): 



~ ~ ~ - _ de 

ftj = psin/3 ^ pfP, (5 = arcsin(Z/p), Z = Z + |, Gj = — 

(JjL) 



9r (15) 

/ 



We have taken into account here the fact that it is the small (3 that are significant for the 
oscillations. In this limit, we assume a power-law dependence of the coupling constant on (5: 

93 = gP'-'- (16) 

In terms of the variables p and P the distribution (13) takes the form 

{1, p-jf < -^gpfP'', 

^-fS^' -ypfP'<p'pf<igpf(^^ (17) 
0, P-Pf> ypfP''- 

Thus, in terms of the distribution function w{p,l), the model (9)-(ll) turns out to 
be the simplest of all the conceivable ones. It represents, in fact, a linear interpolation 
between regions where the occupation numbers of the quasiparticles assume the limiting 
values zero and unity, see Fig. 4. Instead of analyzing the expression for the second variation, 
it is easier to verify the stability of the distribution (13), (17) directly, by effecting all 
the conceivable transfers of the pair of interacting quasiparticles between the states, and 
calculating the corresponding change of the energy of the nucleus. According to Fig. 4, four 
types of quasiparticle transfer are permissible: I-II, I-III, II-II, and II-III. In all cases we 
have for the change of the system energy 

AE > Gf + Gf (18) 

i.e., the distribution of the quasiparticles over the states is stable. To calculate the oscillating 
part of the nuclear energy, we use the result obtained in [1]: 

= ^ (£;- + £;-*), e'''''''^''''+^^E{n,l)dndl, (19) 

where ^ ^ 

E{n, I) = 4:eflw = -4:epw sin p ^ -Aepfpw (20) 

is the energy of the quasiparticles having the corresponding values of the quantum numbers; 
the additional doubling is due to the spin of the individual quasiparticles, to which an energy 



*Thc latter is confirmed by an analysis of the dependence of the spin-orbit interaction on nuclear dimen- 
sions (see, for example, [2]), i.e., in final analysis, on p/ = kfR. It must be stipulated, however, that this 
pertains for the time being to the calculation of the effects in the pf scale. Subsequently, the spin-orbit 
interaction is actually taken into account in a most substantial manner when the final transition is made to 
the physical A^-scale; see [1]. 
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value Ef = —e is assigned (we shall discuss the validity of the last assumption at the end of 
the section). According to [1] (see also [4]) 



27r(2n + /) ^ -in + 4p + 2p(3^, dndl = -pdp cos^ (3d(5 = ^dpd(3. (21) 

TT TT 

Taking also (17) into account, substituting in (19), and performing the simple integration 
with respect to p, we obtain 

= f eM^imP') {eM^i^gPfP') - eM-^i'^gPfP')) (22) 

To determine which values of the exponent k are compatible with the character of the 
experimental data, it is easiest to turn to the limiting case when the convergence of the 
integral (22) is due principally to that term of the argument of the exponential which is 
proportional to g ^. We expand exp(2ii/p//?^) in a series, confine ourselves to the first two 
terms, and substitute in (19); this yields 

^ epf \ {2gpff-*l^ sinydy ^ cos Ai^pf 

' ^kg\ g io y'-'^' 4t 

CO ^ (23) 

cosydy \- sm Aupf I k/2-1 



The singularities of (23) are governed by trigonometric series that can be easily investi- 
gated. A singularity of the type of a finite jump in the derivative dEi/dpf (see Figs. 2 and 
lb) can be obtained only from an even scries in t = 4p/ — 27ip with cosines. When k < 2 
we have 4/k > 2, and the scries of the derivatives diverges uniformly in accordance with 
the Weierstrass criterion. Thus, there are no singularities of this type if /c < 2. ^ To the 
contrary, if /c > 2 the sum of the series exhibits stronger singularities. To verify this, we 
consider the derivative of the series of interest to us; this derivative is given by 



00 . 
sm ut 



/ , a < 1. 



^Strictly speaking, however, the coupling constant can likewise not be exceedingly large, for then the 
situation reduces to effects due to individual nucleons at the s or p levels; see Figs. 3 and 4. This would 
actually be manifest by a complete disappearance of the magic numbers, and in any case the macroscopic 
theory of shell and magic phenomena developed in [1] and in the present paper would no longer hold. As a 
result we obtain the requirement g <C 

^The case k = calls for a special analysis. Although the expression under the summation sign becomes 
somewhat more complicated in this case, the absence of a jump in the first derivative of E, can be easily 
proved (see Appendix 2). We note that the choice A; = is, in a certain sense, of special interest, since it 
corresponds to a Fermi-distribution smearing band whose width is independent of the orbital momentum (this 
case was already mentioned in the introduction; see footnote 3). For example, in the Cooper phenomenon, 
the angular momentum I of the quasiparticle with respect to the geometric center would not play any special 
role, and only the angular momenta of the components of the Cooper pair with respect to their common 
center of gravity would be of importance. We defer additional comments to the last section. 
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Near the singular point t = we have sini/i ~ vt up to a certain hmiting value u — D l/\t\. 
Replacing the summation by integration, we obtain the estimate 



sm i/r 



i/=i 



y ^ ^ W ^^-du ti>^- ^ (24) 

i/« ./o 



Consequently, at /c > 2 we have q; = 4/A; — 1 < 1 and the derivative of the first term of (23) 
tends to infinity as t — >• (the result (24) agrees with rigorous mathematical theorems; see, 
for example, [5]). Thus, only 

A; = 2 (25) 

is compatible with the experimental data. 

To get rid of the lower bound g ^ 1 on the coupling constant, we substitute (25) in (22) 
and integrate; as a result we obtain ultimately 

E,{N, Z) = -^l {/i(^)97l(p;) + UgmPf)} , (26) 



where 



l+g 



/2(^) = ^%-l), (27) 
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9{x) = 1 at a; > and 9{x) = at a; < 0. Plots of the functions /i and /2 are given in Fig. 5. 
The trigonometric series 

^(Pf) = z2 = Y-27r Pf-^P +A(^pf--p) , -{p-l)<pf<-(p + l) 

u=l 

u=l 

describe shell oscillations in spherical nuclei in the presence of a residual interaction ^. Dif- 
ferentiating (26), we obtain the discontinuity 

^s^e^hig) (29) 

of the nucleon binding energy in the nucleus. According to (27), this discontinuity increases 
without limit as g ^ 1. The physical qualitative picture was already explained in the 
introduction and illustrated in Fig. 3. Now, however, unlike in (6), the oscillations tend to 
infinity much less rapidly, namely, logarithmically. The weakening of the singularity is due to 
the presence of region II, which is not free of fermions (see Fig. 4); on the boundary between 
regions I and II, the occupation numbers (17) do not change jumpwise in the considered 
model. 



^Greatest interest, however, attaches nevertheless to the nearest vicinity of the magic nucleus, in which 
the function 01 tends to zero like fln|f|. At the present time it is still not quite clear whether allowance 
for the small term proportional to 01 is an exaggeration of the macroscopic accuracy with which the entire 
developed concept holds (see also [1]). This has no effect whatever on the subsequent results pertaining to 
the magic jump As of the nucleon binding energy. 
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It is quite doubtful whether such a sharply delineated boundary of the region of inter- 
mediate values of occupation numbers actually exists, and actually the coefficient of the 
function D)l{pf) should not become infinite anywhere. We note in this connection the curi- 
ous possibility of constructing a simple interpolation formula for the function fi{g). Since 
the transition gi — > to the case (4) where there is no residual interaction yields /i — > 1, 
and since in the asymptotic region g ^ 1 we have according to (27) /i = i/g^, the following 
interpolation comes to mind: 

Mg)^ 1/(1 + g'). (30) 



We then have, in particular. 



O;^ = (31) 



1 



Thus, the characteristic u; of the deviation from the ideal case (2) and (4), which was 
introduced earlier in [1] from intuitive and empirical considerations, turns out to be di- 
rectly expressed in terms of the constant of the residual interaction between the quasiparti- 
cles. Nevertheless, it seems more natural to start from the very beginning from analytic or 
nearly-analytic expressions for the distribution function w{p,P). One such example will be 
considered in the next section. 

The prospects for choosing the most successful interpolation may become clearer if we 
turn finally to the assumption (made above in connection with formula (20)) that each quasi- 
particle has a definite energy —e. There exists also the energy (14) of the interaction between 
the quasiparticles. When the spectral component E^^^^ of the energy of this interaction is 
calculated by means of formula (19), integration over the region I will certainly not yield 
a macroscopic contribution. Indeed (see formula (14) and the text immediately following 
it) the only inaccuracy in the determination of bj corresponds precisely to the error in the 
interaction energy ~ GjQj. We therefore confine ourselves to integration over the interaction 
region II, and simple calculations lead to 



~ dp 



Since de/dp ~ £/Pf, the contribution (32) to the energy of the nucleus is not macroscopic; 
according to the criteria analyzed in [1], it need not be taken into account. 

This circumstance is apparently not accidental, and is by no way a distinguishing feature 
of just this model. From a more general and physically more lucid point of view, this question 
was already analyzed in the introduction above. 



3 Case of analytic distribution function 

The fact that the contribution of the interaction energy to the oscillations is of no importance 
(at a specified quasiparticlc distribution function) enables us to reformulate the problem in 
the manner outlined in the introduction, namely, the oscillating part of the ground-state 
energy of a spherical nucleus is defined by the function w (which is the diagonal part of 
the density matrix relative to the basis of the quasiparticles in question). The form of the 
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function w{p, /?) is by far not as arbitrary as it might seem at first glance. First, it must tend 
rapidly to zero and unity in the asymptotic regions; it is not very likely that this function 
is non-monotonic. Further, the analysis in the preceding section (see in particular (25) 
and (26)) suggests that only a quadratic dependence of the width of the Fermi-distribution 
smearing band on the orbital angular momentum can be reconciled with the experimentally 
observed character of the magic phenomena at more or less arbitrary intensities of the residual 
interaction. Indeed, it is easy to show that for all distributions of the type w{{p — Pf )/l3'^) 
the oscillating part of the energy Ei always reduces to a linear combination of the expression 
9Jt(p/) and ^(p/). It is therefore natural to attempt to use a function of the type of the 
ordinary Fermi distribution (7), but with a modulus that depends quadratically on the angle 

We calculate the spectral component 

= ^p^fe^'"""^ dppex.p{2iupfp^} j w{p, p) exp{4ii/(p - pf))dp (34) 
of the number of states filled with quasiparticles, and in accord with 

oo 

TVi = ^(iV'^ + N""*) (35) 



we use formula (3) (taking into account the conclusions drawn above concerning the role of 
the energy of interaction between quasiparticles, this is equivalent to formulas (19)): 

^■ = -||g{^c(.i-^).e.}. 

We have put here r = {l/4)gpf (the convenience of such a change of notation will become 
clear later on). Separating in the Riemann ( function 



coo ^^z-l^-qu 

iu = 

the real and imaginary parts, we arrive at the trigonometric series (28): 



^ / X in sm(u/ng) , 



2ng Jq sinh(ii/2) cosh^(l/g') 
Plots of the functions ^1(5^) and ^2(5^) are shown in Fig. 6. 
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Just as for the model example (17) and (25) of the preceding section, in the region (7 S> 1 
of the strong residual interaction there predominates, generally speaking, the function 9Jl(pj), 
which enters here with a coefficient 

F2(^)^7r/2^, 5»1. (39) 

With our choice of the normalization of the couphng constant (see above), Eq. (39) coincides 
with the second formula of (27). Continuing this analogy, we note that in the non-analytic 
model (23), (27) there was no term with 01(p/) at all at 5^ < 1. In our case, however, this 
corresponds to exponential smallness of this term if the coupling between the quasiparticles 
is weak: 

F,{g) - g « 1. (40) 

From the more formal point of view, it is curioTis to note that the limiting behavior of (40) 
indicates the presence of an essential singularity at g = 0. One cannot exclude the possibility 
that this mathematical property of the limiting case — > of an absolutely sharp Fermi 
boundary is quite general. 

The magic singularities are due only to the function 9Jl(p/) (see also footnote 7). The 
formulas characterizing these singularities are 

As seen from Fig. 6, the function ^1(5') has a maximum in the region of intermediate values 
of the couphng constant. Thus, the effect of interest to us does not depend monotonically 
on the intensity of the residual interaction^. 



4 Comparison with experiment 

For the 52 magic nuclei considered in the preceding article [1], we obtained the average 
values pertaining to seven experimentally accessible cusps (four neutron and three proton) 
on the mass curve. The results of such a reduction are given in Table I. In accordance with 
the interpolation (38) and (31), ^/uJ is the crudest characteristic of the residual interaction. 
The reason why cu (or y/cu) may turn out not to be an appropriate characteristic of the 
phenomenon of interest to us are given in footnote 8. The coupling constant g' was calculated 
from formulas (27) and (29), and the quantity g corresponding to it in the analytic case (33) 
was determined from relations (38) and (41). 

^This non-monotonicity is an interesting feature of the phenomenon, of which we were unaware during 
the time of writing of the preceding paper [1]; nor is it revealed by rough interpolations of the type (30). 
The increase of the oscillations at intermediate values of g can be qualitatively understood as a sort of 
resonant amplification due to the similarity in the form of the lower boundary of the region of intermediate 
values of occupation numbers and the corresponding Regge trajectory on Fig. 3. In the non-analytic model 
examples (6) and (29) (see 27)) this resonance was extremely sharply pronounced and caused the effect to 
become infinite as g 1. It can be assumed that in the analytic case, too, the presence of a maximum is not 
accidental, and that if the Fermi boundary is not absolutely sharp, g = 0, the oscillations actually reach their 
maximum swing. For the considered distribution (33), this extremum is given hy g = 0.55, Fi^max = 1.20, 
and LOmin = -0.17. 
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de 
dp 



^P-^. — o- (44) 



With respect to the shell oscillations the dependence of the width of the Fermi-distribution 
smearing band on j3 extends only to the angles 

2p/^' ~ 1, (42) 

after which the integrals describing them (see (22), (34)), in any case, converge rapidly. 
According to (33), the corresponding characteristic width in the p scale is 

Tp^T^^ = \gpf^' = y. (43) 

We change over to the energy scale: 

f 

Here R is the radius of the nucleus. We now assume 

1.2 • 10-^^A^/^[cm] (45) 

and, as an estimate, set the effective mass m* of the quasiparticles equal to the mass of the 
free nucleon. The widths Se in MeV are given in Table I. 

For the concrete values of the orbital angular momentum, the width Se of the smearing 
band turns out to be different, since it depends on / = / + 1/2 quadratically. Taking the 
estimate (42) into account, we confine ourselves to not too large values of /, which are 
characteristic of the considered shells. The results are given in Table II. 

All the characteristics of the residual attraction listed in the table, regardless of the 
degree of their accuracy and of the choice of scale, indicate consistently that its intensity 
decreases rapidly with increasing nuclear dimensions ^. 



5 Conclusion and discussion 

The most interesting result, in our opinion, is the dependence of the residual interaction on 
the orbital angular momentum of the quasiparticle: 

ocl' ~ (/ + 1/2)2. ^4g^ 

^The method of coping with the even-odd mass oscillations used in the reduction of the experimental 
data was already described earlier (see [1], footnote 12). We note in this connection that the change in the 
number of protons by 2 is a reflection of the Coulomb energy of the nucleus is reflected in the CoTilomb 
energy of the nucleus and in its derivative with respect to the number of charged particles. Since in fact 
an energy of pure electrostatic origin has no singularity whatever, it is advisable to subtract the "fictitious" 
effect produced by it from the observed running discontinuities Ae. In connection with the introduction of 
such a correction for the Coulomb interaction, the values of u used in the present paper for the proton magic 
nuclei (see Table I) turn out to be somewhat higher than those previously published (see Table II and Fig. 4 
of [1]). This correction affects little the conclusions concerning the course and magnitude of the residual 
interaction of the protons. 
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Table I 



p 


AT" r7 

i\, z 


Neutron mag 


1 * 

ic nuclei 


Proton 


magic nuclei 






\^ 


9' 


9 


6e (MeV) 




9' 


9 


fc (MeV) 


4 


28 


2.4 


2.7 


3.3 


5.2 


1.8 


2.2 


2.6 


3.8 


5 


50 


1.6 


2.0 


2.3 


3.3 


1.5 


1.9 


2.2 


2.6 


6 


82 


1.4 


1.8 


2.0 


2.6 


1.0 


1.5 


1.6 


1.6 


7 


126 


1.1 


1,6 


1.7 


1.9 











Table II (The values of Se are given in MeV.) 



I 


P = 


= 4 


P = 


= 5 


P = 


= 6 


P 


= 7 




5eN 


Sez 


Sen 


Sez 




Ssz 


Sen 


Sez 





0.4 


0.3 






0.14 


0.08 






1 






1.9 


1.5 






0.8 




2 


10.2 


7.5 






3.5 


2.1 






3 






10.2 


8.1 






4.3 





On the one hand, this result is quite natural. Since we are dealing with a scalar effect, it 
should be expressed precisely through the scalar square of the angular-momentum vector. 

According to the estimate (42), the angular-momentum values that play an important 
role are relatively small. From this point of view, (46) can be regarded as the first term 
of a series expansion in powers of the ratio l/pf. Why does the expansion have no zeroth 
term independent of 11 Does this mean that in a nucleus of finite dimensions the residual 
interaction is due only to the additional integrals of the quasiparticle motion, other than 
the energy? Unfortunately, we see no possibility of giving a sufficiently clear-cut and defi- 
nite answer to such questions. Nonetheless, we wish to point out the difficulties that would 
probably be encountered in attempts to reconcile the obtained picture with the assumption 
concerning the Cooper phenomenon in nuclear matter. It is precisely such a phenomenon 
which is characterized by a constant width of the transition region of the statistical distri- 
bution, having the same value for all the quasiparticles located near the Fermi boundary 
However, the presence of a constant component of the width of this band would lead to the 
absence of the experimentally observed magic jumps of the binding energy of the nucleons 
(see Appendix 2, and also footnotes 3 and 6). The available experimental data therefore 
agree more readily with the simpler and natural hypothesis concerning the character of the 
energy spectrum of unbounded nuclear matter as such. It seems likely that it can represent 
a "normal" Fermi liquid with an absolutely sharp Fermi boundary for quasiparticles [7]. 

On the other hand, spherical nuclei of finite radius have a residual interaction of the 
type characterized by (46). One must not forget, however, that it was actually estabhshed 
only on the basis of data that pertained throughout to the nearest vicinity of some magic 
nucleus. The experimental data give the impression that the residual interaction of the 
considered type decreases gradually with increasing distance from the magic nucleus, and 

good illustration is the well-known problem of an almost-ideal Fermi gas with weak attraction. Its 
solution is frequently given in the momentum representation (see, for example, [6]). It is difficult to visualize 

how the purely formal operation of the transition to the ^-representation could change the width of the 
transition band of the Fermi distribution, and furthermore make it dependent on I (see also footnote 6). 
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becomes somehow restructured in final analysis. The thermodynamics of the transition was 
considered in [8]. The fact that the phase transition due to weakening and restructuring of 
the residual interaction causes also the nuclear shape to become non-spherical can hardly 
be regarded as surprising. It was already shown earlier (see [4]) that in the simplest scheme 
without interaction the sphere is absolutely unstable in general for any number of particles. 

We thank 1.1. Gurevich, L.P. Kudrin, I.M. Pavlichenkov, G.A. Pik-Pichak, V.P. Smilga, 
and K. A. Ter-Martirosyan for a discussion of the results. 



Appendix 1 

In the preceding paper [1] (see also [4]) we considered for concreteness a scheme with a wall 
that is impermeable to the quasiparticles and is located at a distance R from the center of 
the nucleus. The roots of the wave equation of the free motion of the particle in the spherical 
region, shown in Fig. 3 above, corresponded to such a boundary condition. We shall show 
that this does not limit the generality of the results pertaining to shell oscillations of the 
volume energy of the spherical nucleus. 

The Bohr-Sommerfeld quantization rule [2] used to determine the eigenvalues is 

f ki{r)dr = 7i{n + 'j). (Al.l) 

J a 

Here 



ki{r) = ^Jk^ -P/r^, 1 = 1 + 1/2, (A1.2) 

and the lower limit of integration is due to the centrifugal barrier; since the nuclear matter 
is homogeneous, the wave number k is constant in the interval region. The additional phase 
7 depends on the properties of the true structure of the transition layer on the surface of 
the nucleus. We change over to the dimensionless variable kr = p': 

J'n{p')dp = 7r{n + ^), n=sJl-P/p'\ (A1.3) 

The Regge trajectories responsible for the oscillations are characterized by the relation 

2n + l=p, p = 2,3,4,..., (A1.4) 
between the quantum numbers (see the Introduction and Fig. 3). Therefore 

dn/dl = -1/2, d^n/dt^ = 0. (A1.5) 
We now differentiate the entire relation ()A1.3|) along the trajectory 



Elementary integration yields 



TX ^ TT d-f 

/cipl— — arccos - = hvr— . (A1.7) 

^'^^dl p 2 dl ^ ^ 
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It follows therefore that the derivative dp/dl vanishes at / = (the last term in the right- 
hand side, which is due to the differentiation of the phase increment, does not affect the 
result (see fjA1.9|) below). In the second differentiation of the Bohr-Sommerfeld formula, we 
omit the terms that are known to vanish at the extremum: 

We now take into account the fact that the phase 7 is determined by the actual details 
of the nuclear interactions in the nearest vicinity of the upper limit of integration. Any 
characteristic that influence the results can depend here on the angular momentum I only 
in the combination l^. Therefore 

S = 2/i, ^ = 2^^\. (A1.9) 
dl dP dP dP p2 

Consequently, at the maximum point / = we have 

d^p/dP = -1/p, (ALIO) 

neglecting the terms ~ Vp^- Thus, the form of Eq. (1) 

Ap = -(/ + l/2)V2p (Al.ll) 

does not depend on the details of the nuclear surface-layer structure. We determined finally 
the ordinates of the extrema of the successive trajectories. AtZ = Z + l/2 = 0we have 
according to ()A1.4jl n = p/2 + 1/4. Substituting in ()A1.4|1 . we get 

Pmax = vr(ip + 7'), (A1.12) 

where 7' = 7 + //4. Replacing pmax by kjR = pf, we see that ()A1.12|1 coincides in fact with 
the previously published (see [1], formula (19), and footnote 9) "quantization rule" for the 
magic values of this parameter. 



Appendix 2 

We consider here the special case k = 0, for which there is no asymptotic region described by 
formulas of the type (23). The width of the band II of the smearing of the Fermi distribution 
(see Fig. 4), which does not depend on /, will be denoted by go. Formally, this requires that 
we put k = in (17) and make the substitution gpf qq. Substituting then in (34) and 
(35), we obtain after a simple integration 

iV, = ^V^^"'"^°r''^^- (A2.1) 

In the next application of formula (3), we transform also the trigonometric expressions 
under the summation sign 
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Here 

i/=i ■' 

Since the expression for 9?l(p/) (see formula (28) above, and also [1]) is characterized by 
a jump of the first derivative, only the second derivative of the function 2{pf) becomes 
discontinuous. Moreover, these singularities (which are located in this case at the inflection 
points of the function Ei{pf)) are now shifted by a distance ±go/2. The physical reason for 
this is seen directly from Fig. 4. At the very point on the pf scale where the first derivative 
would experience a discontinuity at = 2, the expression ()A2.2jl has an ordinary analytic 
minimum. Thus, the nucleon binding energy e experiences no discontinuity anywhere. 
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Figures captions 

Fig. 1. 

Fig. 2. Schematic diagram of the shell oscillations of the nuclear mass in the vicinity of 
the magic number N = 126. 

Fig. 3. Roots of the Bessel functions with half-integer index on the (/,p) diagram. 

Fig. 4. 

Fig. 5. 

Fig. 6. 
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